ESC-TR-2001-074 


Project  Report 
FPR-14 


Analysis  of  SAR  Image  Formation  Equations 
for  Stationary  and  Moving  Targets 


A.F.  Yegulalp 


20  June  2002 


Lincoln  Laboratory 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
Lexington,  Massachusetts 


Prepared  for  the  Defense  Advanced  Research  Projects  Agency 
under  Air  Force  Contract  F19628-00-C-0002, 

Approved  for  puhhc  release;  distribution  is  unlimited. 


20020705  059 


This  report  is  based  on  studies  performed  at  Lincoln  Laboratory,  a  center  for 
research  operated  by  Massachusetts  Institute  of  Technology.  This  work  was  spon*- 
sored  by  the  Defense  Advanced  Research  Projects  Agency,  IXO,  under  Air  Force 
Contract  F1962R>00“C-0002. 

This  report  may  be  reproduced  to  satisfy  needs  of  U.S.  Government  agencies. 


The  ESC  Public  Affairs  Office  has  reviewed  this  report,  and 
it  is  releasable  to  the  National  Technical  Information  Service,  where 
it  will  be  available  to  the  general  public,  including 
foreign  nationals. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


FOR  THE  COMMANDER 


Contracted  Support  Management 


Non-Uncoln  Recipients 

PLEASE  DO  NOT  RETURN 
Permission  is  given  to  destroy  this  document 
when  it  is  no  longer  needed. 


Massachusetts  Institute  of  Technology 
Lincoln  Laboratory 


Analysis  of  SAR  Image  Formation  Equations  for 
Stationary  and  Moving  Tai^ets 


AS.  Yegulalp 
Group  45 


Project  Report  FPR-14 


20  JuBe  2002 


Approved  for  public  release;  distribution  is  unlimited. 


Lexington 


Massachusetts 


ABSTRACT 


This  report  presents  a  detailed  mathematical  analysis  of  backprojection  and  “o)-k” 
(range  migration)  image  formation  equations  for  synthetic  aperture  radar  (SAR).  The 
image  processor  transfer  fimction  is  computed  for  both  stationary  and  moving  targets, 
including  antenna  and  target  pattern  effects.  Assuming  straight-line  platform  motion,  it  is 
demonstrated  that  both  image  formation  algorithms  provide  exact  focusing  of  stationary 
targets.  Moving  targets  are  displaced  to  false  locations  and  blurred  into  elliptic  or 
hyperbolic  curves.  Algorithms  are  suggested  which  will  exactly  refocus  moving  targets 
in  SAR  images  based  on  only  one  unknown  motion  parauKter.  All  results  are  derived 
without  narrow-band,  narrow-angle,  or  narrow-swath  assumptions. 

Sinular  derivations  for  moving  target  signatures  and  focusing  algorithms  have  been 
presented  in  [1].  Reference  [1]  also  demonstrates  the  application  of  moving  target 
focusing  to  experimental  data. 
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1.  INTRODUCTION 


1.1  SAR  GEOMETRY 

Figure  1  shows  the  basic  geometry  for  a  synthetic  aperture  radar  (SAR)  data  collection. 
The  radar  platform  is  assumed  to  follow  a  straight-line  trajectory,  with  its  antenna  beam 
illuminating  terrain  to  one  side  of  the  aircraft.  Points  on  the  terrain  surface  may  be  labeled 
according  to  the  along  track  (cross-range)  coordinate  x  and  the  cross-track  (slant-range) 
coordinate  y.  The  so-called  “slant-plane”  coordinates  {x,y)  are  natural  for  image  formation,  and 
will  be  used  throughout  this  report.  Note  that  y  is  by  definition  positive.  Converting  from  the 
slant-plane  coordinates  to  fixed  “real-world”  coordinates  (e.g.,  latitude-longitude-elevation)  is 
straightforward,  but  does  require  precise  knowledge  of  the  terrain  contour  and  the  flight  track 
geometry. 


Figure  1.  SAR  flight  geometry. 
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In  practice,  the  aircraft  will  deviate  somewhat  from  a  straight-line  path.  Both  image 
formation  algorithms  studied  in  this  report  can  be  adapted  to  compensate  for  off-traek  motion, 
but  only  the  baekprojection  algorithm  can  provide  exact  image  focusing  for  arbitrary  aircraft 
motion  and  ground  topography.  To  reduce  the  scope  and  complexity  of  the  analysis  in  this 
report,  all  calculations  will  be  performed  for  strietly  straight-line  motion. 

Figure  2  shows  the  radar  platform  moving  past  a  stationary  target  on  the  ground.  A  fixed 
antenna  beam  is  assumed,  illuminating  a  region  of  angular  width  6ant  on  the  ground.  Most  of  this 
doeument  will  assume  that  the  beam  is  uniform  within  the  angle  6ant  and  zero  outside.  Section  4 
will  show  how  arbitrary  antenna  patterns  may  be  accommodated  using  a  simple  modification  of 
the  uniform  beam  results. 


Full  Synthetic  Aperture  for  Entire  Image 
Aperture  Relevant  to  Target 


Aperture  Processed  for  Target 


Figure  2.  Beam  and  aperture  geometry. 
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From  the  point  of  view  of  the  target,  the  platform  sweeps  out  an  aperture  of  angular  width 
dant  during  the  time  it  illmninates  the  target.  The  SAR  processor  may  use  any  part  of  the  full 
aperture  in  forming  an  image. 

In  stripmap-mode  SAR,  the  processed  aperture  corresponds  to  a  user-selected  integration 
angle  Omt  ^  dant-  Cross-range  resolution  in  the  image  improves  as  6^  is  increased.  Using  a 
constant  integration  angle  throughout  the  image  produces  an  image  with  spatially  invariant 
resolution  and  sidelobe  structure.  Note  that  the  processed  aperture  is  referenced  with  respect  to 
the  point  being  imaged,  so  different  points  in  the  image  will  be  formed  using  different  portions 
of  the  full  aperture.  Moreover,  the  length  of  the  processed  aperture  will  be  proportional  to  range 
in  order  to  maintain  a  constant  integration  angle. 

Spotlight-mode  SAR  uses  all  available  aperture  data  to  form  the  image.  Spotlight-mode 
radars  also  may  employ  a  steerable  beam  to  maintain  illumination  on  a  fixed  patch  on  the  groimd 
over  a  longer  aperture.  For  clarity  of  exposition,  this  document  will  focus  on  stripmap-mode 
SAR,  but  most  results  (aside  firom  antenna  pattern  effects)  extend  to  spotlight-mode  in  a 
straightforward  manner. 

1.2  SAR  DATA  MODEL 

1.2.1  “Hopping”  Model 

As  the  radar  platform  moves  along  the  synthetic  aperture,  it  transmits  a  series  of  pulses  and 
records  the  resulting  backscatter.  To  simplify  analysis,  a  “hopping”  model  will  be  used  in  which 
the  platform  hops  forward  from  one  pulse  to  the  next,  but  remains  stationary  during  the 
transmission  and  reception  of  each  pulse.  The  approximation  causes  negligible  errors  as  long  as 


7®  <0.1 


c 

\dRldt\’ 


(1) 


where  T  is  the  pulse  lengtfi,  B  is  the  transmit  bandwidth,  c  is  the  speed  of  light,  and  dRjdt  is  the 
range  rate  of  change  for  groimd  points  in  the  radar  beam  [2,3].  Typically  the  constraint  is  easily 
satisfied  for  airborne  platforms,  but  may  be  violated  for  orbital  platforms. 

The  hopping  model  may  also  be  applied  to  moving  targets.  The  target  is  assumed  to  move 
between  pulses,  but  remain  stationary  during  each  pulse.  Equation  (1)  determines  the  validity  of 
file  hopping  model  for  moving  targets,  provided  that  the  range  rate  dRjdt  includes  the  combined 
effect  of  platform  and  target  motion. 
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1.2.2  Range  Compressed  Pulse  Data 


No  particular  assumptions  will  be  made  regarding  the  radar  waveform.  Instead,  it  will 
simply  be  assumed  that  the  return  pulses  are  range  compressed  so  that  an  ideal  point  scatterer’ 
with  unit  cross  section  at  range  R  will  )deld  the  signal 

=  smc{B{t-2R/c)) ,  (2) 


where  B  is  the  radar  bandwidth,  fc  is  the  center  frequency,  and  sinc(«)  =  .  Expressed  in  the 

frequency  domain.  Equation  (2)  is  an  integral  over  a  flat  spectrum: 


1  ymax 

S{t,R)=-^  J  (3) 

^  /min 

The  minimum  and  maximiim  frequencies  fmm  and  fmax  are  related  to  the  radar  bandwidth  and 
center  frequency  by  f^=fc-  Bjl  and  Bjl . 

The  range-compressed  pulse  data  will  be  denoted  by  D(^t),  where  ^  is  the  radar  position 
along  the  aperture  and  t  is  elapsed  time  after  pulse  transmission.  If  there  is  just  a  single 
stationary  point  scatterer  at  (xo,yo)  with  unit  scattering  amplitude,  the  target  range  is 

^(4)  =  yl(xQ-^^  +  yl ,  (4a) 

and  the  range-compressed  data  is 

=  (4b) 

If  the  terrain  is  described  by  a  complex  scattering  amplitude  density  G(x,y),  the  range 
compressed  data  will  be 


D(^,t)  =  l\G(x,y)s{t,^(x-^f  +/ )  dxdy.  (5) 

The  case  of  a  single  point  scatterer  is  equivalent  to  G(x,y)  =  ^(x-Xo)<J(y-yo). 


1 


By  (definition,  an  ideal  point  scatterer  has  constant  scattering  amplitude  as  a  function  of  frequency  and  direction. 
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Equations  (2)  through  (5)  implicitly  assume  uniform  gain  in  frequency  and  direction  for 
both  the  radar  antenna  and  the  terrain  scattering.  The  more  realistic  case  of  nonuniform  gain  will 
be  addressed  in  Section  4. 

1.23  Moving  Target  Model 

Consider  a  point  scatterer  target  moving  at  constant  velocity  (with  respect  to  slant  plane 
coordinates).  Working  within  the  assumptions  of  the  hopping  model  of  Section  1.2.1,  the 
corresponding  range-compressed  data  can  be  written  as 

D„(#,0  =  S(/A(#)),  (6a) 

where  the  target  range  is 

and  the  target  coordinates  are 


The  target  motion  is  parameterized  using  the  broadside  position  (xo,yo)  and  the  normalized 
velocity  components  (ux,Uy).  The  broadside  position  is  simply  the  position  of  the  target  at  the 
moment  when  it  is  directly  broadside  to  the  aircraft  (i.e.,  when  ^x„).  The  normalized  velocity 
components  are  the  target  velocity  components  (VjcVyj  divided  by  the  platform  speed  Vpiatform- 


^platform 


Uy=  - ^ 

^platform 


Straightforward  algebraic  manipulation  shows  that 


(8) 


(9) 


where 


r=^(l-u^f+uj  , 


and 


1-Mv 


(10) 
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The  new  quantities  defined  in  Equation  (10)  all  play  important  roles  in  the  upcoming  results.  It 
will  be  shown  that  the  motion  parameter  ^completely  determines  point  target  defocusing:  at  a 
given  range,  targets  with  the  same  yhave  the  same  blurring.  (Note  that  ^'may  be  understood  in 
physical  terms  as  the  target  speed  measured  in  the  rest  frame  of  the  aircraft  and  normalized  by 
the  aircraft  speed).  It  will  also  be  shown  that  the  quantities  Xs  and  ys  are  the  shifted  (apparent) 
coordinates  of  the  target  in  the  image. 

For  later  notational  convenience,  the  parameter  e  is  defined: 


1  _  — 2w^+i/^ 
\-2u^+u^ 


(11) 


Stationary  targets  correspond  to  y  =  1  and  e  =  0.  Both  y  and  e  will  be  used  freely  in  writing 
results,  but  it  should  be  remembered  that  they  are  not  independent  quantities. 


13  SARAS  A  LINEAR  SYSTEM 


The  process  of  SAR  data  collection  can  be  thought  of  as  a  linear  mapping  Ldc  that 
transforms  the  terrain  scattering  amplitude  into  the  range  compressed  data: 

The  transformation  is  linear  in  the  sense  that  a  linear  combination  of  scattering  amplitudes  will 
map  into  a  linear  combination  of  range-compressed  data: 

LDc{flr<j,(j:,y)  +  yff(72(x,y)}  =  flrLDc{Gi(Jc,y)}  + A^dc{<J2(^.7)}  • 


SAR  image  formation  can  also  be  thought  of  as  a  linear  transformation  Lip  .  It  maps  the 
range-compressed  data  into  the  image  data,  attempting  to  recover  the  original  terrain  scattering 
amplitude: 


I(x,y)  =  L„  {23(1.0}  =L,f  {G(x,:f)}}  =  G(x,y) . 


The  recovery  can  only  be  approximate  due  to  the  limitations  of  finite  radar  bandwidth  and 
aperture.  In  general,  the  image  is  a  filtered  version  of  G(x,y).  The  properties  of  any  given  image 
formation  algorithm  can  be  completely  characterized  by  computing  the  transfer  function  of  the 
aggregate  transformation  LifLdc.  By  linearity,  it  suffices  to  compute  the  image  function 
resulting  from  a  point  scatterer  with  unit  scattering  amplitude. 


The  analysis  presented  in  Sections  2  and  3  will  focus  on 
spectrum  i{k^  ,ky)= \\dxdy  I (x,  .  Starting  from  the 


computing  the  image  spatial 
stationary  point  scatterer  data 
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in  Equation  (4),  it  will  be  shown  that  both  the  backprojection  and  co-k  image  fonnation  equations 
produce 


IiK,ky)  = 


0  fox  (k^,ky)i.  Q 


The  band-support  region  Q.  is  illustrated  in  Figure  3,  and  defined  as  follows: 

Q  =  \  {k^,ky)  I  <^k<f^,  > 0  and  y  < sin%- 


where  k 


^+kl  . 


(12) 


(13) 


ky 


Figure  3.  Band  support  region  £2 for  SAR  image. 


Equation  (12)  shows  that  the  backprojection  and  co-k  image  formation  equations  have  ideal 
bandpass  transfer  functions.  In  other  words,  the  image  is  a  bandpassed  version  of  the  original 
scattering  amplitude  G(x,y),  with  no  extra  phases  or  amplitude  factors  introduced.  The  size  and 
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shape  of  the  band-support  region  Q  reflects  the  transmitted  waveform  bandwidth,  carrier 
frequency,  and  processed  integration  angle. 

1.4  MOVING  TARGETS  IN  SAR  IMAGES 

One  of  the  basic  assumptions  in  standard  SAR  image  formation  is  that  targets  and  clutter 
do  not  move  while  being  illuminated  by  the  radar  beam.  If  a  target  does  move,  it  will  generally 
appear  blurred  and  displaced  to  a  false  location  in  the  image.  The  analysis  in  Sections  2  and  3 
will  consider  a  moving  point  target  and  derive  the  resulting  image  in  both  the  spectral  and  spatial 
domains.  In  the  spectral  domain,  the  moving  target  image  data  looks  like  the  stationary  target 
result  of  Equation  (12)  multiplied  by  a  nonlinear  phase  function  and  a  nonconstant  amplitude 
function.  In  the  spatial  domain,  the  moving  target  is  blurred  into  an  elliptic  or  hyperbolic  arc. 

The  derived  moving  target  signature  may  be  used  as  the  basis  of  an  algorithm  that 
reprocesses  a  standard  SAR  image  to  focus  moving  targets  (and  consequently  blurs  stationary 
targets).  Despite  the  fact  that  the  target  velocity  has  two  independent  components,  it  will  be 
shown  that  there  is  only  one  effective  motion  parameter  that  controls  moving  target  focusing.  In 
other  words,  an  algorithm  that  seeks  to  focus  moving  targets  of  unknown  velocity  need  only 
search  over  a  o«e-dimensional  space  of  possible  focusing  parameters.  A  second  parameter  is 
required  only  for  precise  compensation  of  motion-related  distortions  in  the  target  and  antenna 
patterns. 

A  closely  related  discussion  of  moving  target  focusing  is  presented  in  [1]. 
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2.  BACKPROJECTION  IMAGE  FORMATION 


2.1  BACKPROJECTION  EQUATION 

Given  range-compressed  data  D(^t),  define  the  baclqjrojection  preprocessed  data  by 


iTti 


Att 


dt 


(14) 


The  backprojection  image  is  given  by  the  equation 

I(x,y)  =  y  i  F„(4,i^(x-if+y^]d^  (15) 

x-ytan^- 

The  basic  baclq)rojection  algorithm  has  been  discussed  in  a  number  of  places  in  the  SAR 
literature  [3,4].  The  exact  details  of  the  baclq)rojection  equation  vary  somewhat  among  different 
authors;  the  preceding  definitions  are  chosen  to  yield  an  ideal  bandpass  transfer  function. 

In  practice.  Equations  (14)  and  (15)  are  applied  to  digitally  sampled  data.  The  integral  in 
Equation  (15)  is  approximated  using  a  sum  over  discrete  aperture  positions.  Computing  the 
integrand  in  Equation  (15)  requires  a  suitable  interpolation  in  t  for 

The  derivative  in  Equation  (14)  must  also  be  adapted  to  discrete  data.  The  easiest  approach 
is  to  combine  the  derivative  with  the  range  compression  process.  Frequently  compression  is 
performed  by  convolving  a  matched  filter  M(t)  with  raw  pulse  data  The  derivative  is 

given  by 

In  other  words,  the  derivative  can  be  applied  just  once  to  the  matched  filter  M(t)  prior  to  any 
processing.  The  range  compressed  data  will  then  automatically  be  differentiated,  with  no  extra 
computation  required.  The  derivative  of  a  discretely  sample  matched  filter  can  be  performed 
easily  in  the  frequency  domain  by  multiplying  by  a  factor  of  frequency. 

In  recent  years,  a  number  of  accelerated  backprojection  algorithms  have  been  developed 
[5-7].  The  accelerated  algorithms  do  not  change  the  basic  structure  of  Equations  (14)  and  (15). 
Instead,  they  break  the  required  numerical  calculations  into  stages  of  increasing  resolution, 
organized  to  reduced  the  total  number  of  operations.  The  analysis  in  this  section  should  apply 
equally  well  to  any  of  the  accelerated  backprojection  algorithms. 
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2.2  STATIONARY  PHASE  ANALYSIS  (STATIONARY  TARGET)^ 


Assume  that  the  radar  illuminates  a  stationary  point  scatterer  target,  as  described  in  Section 
1.2.2.  The  range  compressed  data  is  given  by  Equation  (4b).  Plugging  the  range  compressed 
data  into  Equation  (14)  yields  the  backprojection  preprocessed  data: 


^  \n 


J  f 


df. 


(16) 


Plugging  Equation  (16)  into  the  backprojection  Equation  (15)  and  computing  the  spatial 
spectrum  yields  the  four-dimensional  integral 


hK,k,)Mdy] 

^0  ^ 


/mar 

dx  J  d^  df 

6‘  f 

x-ytzn-^ 


yf  -ixkx-iyky  -Rq  (f 


(17) 


The  remainder  of  this  section  will  be  devoted  to  evaluating  Equation  (17)  using  the  method  of 
stationary  phase^. 

The  first  step  is  to  make  the  following  change  of  variables: 

A-tt 

(0  =  — /,  =  and  C  =  (18) 

c 

The  spatial  spectrum  becomes 

oo  ytan^  ^  -^/max 

hk„k,)=-ye-‘^-fdy  J  dzjdi  I  (19) 

0  f  +>0 

where  the  phase  function  <E>  is  defined  as 

^(z, yX,^) = -#x  -  +  yo  ■  (20) 


The  first  derivatives  of  <E>  are 


2 

The  stationary  target  case  is,  of  course,  subsumed  within  the  moving  target  case  treated  in  Section  2.4.  The  reader 
may  regard  the  stationary  target  calculation  as  a  “warm-up”  for  the  more  general  discussion  to  come. 

^  See  Appendix  for  a  brief  review  of  the  multidimensional  stationary  phase  approximation. 


10 


Solving  the  equation  VO  =  =  0  shows  that  there  is  a  single  stationary  point 

Qr,y,^,6>)  =  ,  where 

f  =  i^  =  JoSgn(^^),  C  =  -yQ^y  and  d)=k  =  y[k[+^ .  (22) 

(Note  that  >^0  >  0  since  slant  range  is  an  intrinsically  positive  quantity).  If  the  stationary  point 

falls  outside  the  integration  limits  of  Equation  (19),  the  integral  is  zero.  A  simple  analysis  shows 
that  the  stationary  point  is  inside  the  limits  of  integration  if  and  only  if  (k^,ky)e  Q ,  where 

fi=j 


The  band  support  region  Q  is  illustrated  in  Figure  3.  The  matrix  of  second  derivatives  of  O  is 


A  tedious  but  straightforward  computation  gives  the  determinant: 


11 


(25) 


det9^a)  =  - 


2  2 

yp^ 


ylz^+y^C^+ypf^ 

Plugging  in  the  values  from  Equation  (22)  for  the  stationary  point  yields 


d^6  =  d^0(Z,y,C,d))  = 


k^ 

IXy 

kK 

0 

ypk^ 

k 

Kk] 

kX 

0 

K 

yok" 

ypk^ 

k 

kl 

k 

0 

0 

y 

ypk^ 

k 

K 

K 

0 

K  k 

k 

k 

y 

and 


det926  =  -^ 

ypk^ 


(26) 


(27) 


The  determinant  is  positive  definite  since  ^0  in  the  band  support  region  Q. 

Consequently,  cannot  have  any  zero  eigenvalues.  A  numerical  eigenanalysis  with 

k^=ky  =  yQ=l  shows  that  there  are  two  positive  and  two  negative  eigenvalues.  Any  continuous 

deformation  of  the  parameters  kx,  ky,  and  yo  will  results  in  a  continuous  deformation  of  the 
eigenvalues.  Since  eigenvalues  cannot  cross  zero,  it  follows  there  will  always  be  two  positive 
and  two  negative  eigenvalues.  Consequently, 


sigd^6  =  2-2  =  0. 


(28) 


Combining  the  information  from  Equations  (22),  (27),  and  (28)  with  the  stationary  phase 
Equation  (A.8)  leads  to  the  result 


iyoky 

0  fox  (k^.ky)€Q 


(29) 
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Equation  (29)  shows  that  the  backprojection  image  formation  equation  has  an  ideal 
bandpass  transfer  function.  In  other  words,  the  image  is  a  bandpassed  version  of  the  original 
scattering  amplitude  G(x,y),  with  no  extra  phases  or  amplitude  factors  introduced.  The  size  and 
shape  of  the  band-support  region  Q  reflects  the  transmitted  waveform  bandwidth,  carrier 
frequency,  and  processed  integration  angle. 

23  IMPACT  OF  PREPROCESSING  STEP 

It  is  possible  to  form  images  without  performing  the  preprocessing  step  in  Equation  (14). 
In  other  words,  the  image  is  computed  directly  from  the  range-compressed  data: 


x+yxan 


^nt 


(30) 


x-ytan- 


^nt 


Equation  (30)  may  be  analyzed  exactly  as  in  Section  2.2.  The  only  change  is  that  the  integrand 
in  Equation  (17)  loses  a  factor  of/ and  some  multiplicative  constants.  Since  the  phase  function 
O  is  the  same,  the  coordinates  of  the  stationary  point  in  Equation  (22)  and  the  matrix  of  second 
derivatives  in  Equation  (24)  remain  unchanged.  The  final  result  is 


for <k„k,)BS 
0  {or(k^,ky)€Q 


(31) 


Comparing  Equations  (31)  to  (29)  shows  that  skipping  the  preprocessing  step  introduces  a 
non-uniform  gain  factor  proportional  to  % .  For  ultra-wideband  radar,  the  extra  non-uniform 
gain  factor  can  have  a  measurable  impact  on  resolution  and  sidelobe  structure  in  the  image. 


2.4  MOVING  TARGET  ANALYSIS 

2.4.1  Stationary  Phase  Calculation  of  Image  Spectrum 

Consider  a  point  scatterer  target  moving  at  constant  velocity,  as  modeled  in  Section  1.2.3. 
The  moving  target  data  may  be  processed  into  an  image  using  the  standard  backprojection 
Equations  (14)  and  (15).  The  resulting  image  will  be  analyzed  in  this  section  using  the  same 
stationary  phase  analysis  applied  to  the  stationary  target  case  in  Section  2.2.  Putting  the  range- 
compressed  data  from  Equation  (6a)  into  the  preprocessing  Equation  (14)  yields 
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(32) 


FBp{^,t)  = 


An 


k 


cK(i) 


7  / 

/min 


Plugging  the  preprocessed  data  from  Equation  (32)  into  the  backprojection  Equation  (15) 
and  computing  the  spatial  spectrum  yields  the  four-dimensional  integral 


Iik,,ky)  =  -j]dy]dx  J  d^\df-^e  ^  ^ 


x-y\ari^-  Aiin  m(^) 


■  (33) 


Note  that  Equation  (33)  is  the  same  as  the  stationary  target  case  in  Equation  (17),  except  for  the 
substitution  The  remainder  of  this  section  will  be  devoted  to  evaluating  Equation  (33) 

using  the  method  of  stationary  phase'*. 

The  first  step  is  to  make  the  following  change  of  variables: 


0)  =—f,  Z  =  x-^,  and  C  =  r(^-x,). 
c 


(34) 


The  spatial  spectrum  becomes 


ytan^ 


An  f 
^  /max 


/  0  Mr  .  b  -^ys 

f.  ^min 


0  -ytan^^ 


(35) 


where  the  phase  function  ^  is  defined  as 


^iZ,y,C,ii>)= -XK  -^-yky+  W/+/  - 0)W~^s  • 


(36) 


The  first  derivatives  of  O  are 


4. 

-X  “ 

dx 


^/7v 


S9__k, 


Co) 


r 


4 


See  Appendix  for  a  brief  review  of  the  multidimensional  stationary  phase  approximation. 
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(37) 


30  _  y(0 


Solving  the  equation  VO  =  =  0  provides  the  stationary  point  coordinates 


z=ys 


y=ys 


(38) 


where  e  =  \  — ^ . 

r 

If  the  stationary  point  falls  outside  the  integration  limits  of  Equation  (35),  the  integral  is 
zero.  A  simple  analysis  shows  that  the  stationary  point  is  inside  the  limits  of  integration  if  and 
only  if  (k^,ky)e  Q y,  where 

(39) 

If  ^>sin(^i„(/2),  the  band-support  region  Qy  is  the  same  as  in  the  stationary  target  case 
illustrated  in  Figure  3.  If  the  angular  width  of  the  band  support  region  is 

reduced  as  shown  in  Figure  4. 

The  matrix  of  second  derivatives  of  <E>  is 


x^^y) 


/min  < < /max .  > 0  and  ^ <  min(sin%i-, y) 
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Figure  4.  Band  support  region  £2y  for  moving  target  SAR  image  when  f  <  sin 
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A  tedious  but  straightforward  computation  gives  the  determinant: 


detd^O  =  - 


yW 


\3/2  * 


Plugging  in  the  values  from  Equation  (38)  for  the  stationary  point  yields 


a"«l>=a'a>(i-,y,C6>)= 


k]^k]+ekl 

yJ^^ 

kxkyyjky+sk^ 

ysk^ 

0 


-X. 

k 


and 


k^kyyjky  +sk^ 

y.k^ 

0 


.  {k]+Eklf 
deta^<I)=  ^ 


yy 


•'X 

k 


0 

ysk^ 

ky 


ky 

0 


(41) 


(42) 


(43) 


It  may  be  easily  shown  that  the  determinant  is  positive  definite  in  the  band  support  region 
iiy.  Consequently,  a^4>  cannot  have  any  zero  eigenvalues.  A  numerical  eigenanalysis  with 
]c^=ky  =  y^=l  and  £■  =  0  shows  that  there  are  two  positive  and  two  negative  eigenvalues.  Any 
continuous  deformation  of  the  parameters  kx,  ky,  y*  and  f  will  results  in  a  continuous  deformation 
of  the  eigenvalues.  Since  eigenvalues  cannot  cross  zero,  it  follows  there  will  always  be  two 
positive  and  two  negative  eigenvalues.  Consequently, 

siga^4>  =  2-2  =  0.  (44) 

Combining  the  information  from  Equations  (38),  (43),  and  (44)  with  the  stationary  phase 
Equation  (A.8)  leads  to  the  result 
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(45) 


l(k^,ky)  =  < 


ky  ^skj[  iy^  ^ky+ £k^ 


Uky-^Sk^ 


r{k 

0 


for  (k^,ky  )E  Qy 
for  (k^,ky)i  Qy 


The  spectrum  depends  on  the  target  motion  purely  through  the  single  motion  parameter 
Y  =  +Uy  .  (Recall  that  the  parameter  fis  a  function  of  y,  as  defined  in  Equation  (1 1)). 

2.4.2  Moving  Target  Signature  in  Spatial  Domain 

The  moving  target  signature  in  the  spatial  domain  is  given  by  the  inverse  Fourier  transform 
of  the  spectrum  in  Equation  (45): 


liX,y)=lldk^dky  hk^,ky)e‘^^*‘^^^  =jjdk^dky  — 


kj+Sk^ 


(46) 


The  phases  have  been  collected  into  the  function  T': 

^ik^,ky)  =  (x-xjk^  +  yky  - yyyjky  +  et 
The  derivatives  of  the  phase  function  are 


(47) 


^  =x-x-y.  , 

9*, 


ek^ 


and  *7 —  ^  y~ys - - 


dk,, 


(48) 


Solving  the  stationary  phase  equation  VT  =  =  0  shows  that  there  are  no  stationary 

X  y 

points  unless 


j(x-Xyf  +  y^  =  yy  and  y>0.  (49) 

When  Equation  (49)  holds,  the  stationary  phase  equation  is  solved  by  all  points  (k^,ky) ,  which 
satisfy 
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(50) 


K  _x-Xs 

ky  ey 

In  other  words,  the  solution  is  a  “stationary  line”  in  the  (kj^,ky)  plane  rather  than  a  single 
stationary  point.  Consequently,  the  stationary  phase  equations  in  the  Appendix  do  not  directly 

apply  to  Equation  (46).  The  existence  of  a  stationary  line  implies  that  det  =  0 ,  so  a  formal 
evaluation  using  Equation  (A.  8)  will  produce  a  meaningless  infinite  answer.  Nevertheless,  the 
basic  principle  underlying  the  stationary  phase  method  still  applies  to  Equation  (46);  the  integral 
should  be  “large”  when  stationary  points  exist  and  small  otherwise.  It  follows  that  the  blurred 
target  will  appear  as  a  curve  in  the  (x,y)  plane  satisfying  Equation  (49).  For  £>0  the  curve  will 
be  elliptic,  and  for  €<0  it  will  be  hyperbolic.  Equation  (50)  determines  the  length  of  the  curve 
since  the  stationary  points  (k^,ky)  are  required  to  lie  in  the  band  support  region  iiy. 

The  curve  of  the  blurred  target  signature  may  be  written  in  parametric  form  as  follows: 

For  £■> 0  (i.e.,  y>l): 

X  =  +  Vf  y,  sin  (p , 

y  =  y^co%(p, 

where  |tan (p\<4e min |sin j .  (51) 


For  £•  < 0  (i.e.,  y<l): 


y-y^  cosh  (p , 

where  |tanh  min  |  sin  ,  y  j .  (52) 

The  curves  in  Equations  (51)  and  (52)  are  illustrated  in  Figure  5.  Note  that  curves  with  £>  0 
always  bend  toward  the  flight  track,  while  curves  with  £<0  curve  away.  The  point  (x^,y,)  is 

the  center  point  of  the  blurred  target  curve;  it  is  related  to  the  true  target  coordinates  via 
Equation  (10). 
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(■*^>y)  ~  e  <  0  (hyperbolic) 


Figure  5.  Moving  target  signature  in  spatial  domain. 


2.4.3  Expansion  for  Small  Target  Velocities 

When  « 1  and  « 1 ,  the  following  leading-order  expansions  become  relevant: 


and 


Fs 


(53) 


At  leading  order,  defocusing  (governed  by  e)  depends  only  on  the  along-track  velocity 
component.  Target  displacement  depends  only  on  the  cross-track  velocity  component,  and  there 
is  no  cross-track  displacement  at  first  order. 

2.4.4  Refocusing  Moving  Target  Images  Via  Spatially  Variant  Filtering 

As  demonstrated  in  the  previous  sections,  a  moving  target  will  appear  as  an  elliptic  or 
hyperbolic  streak  in  a  standard  SAR  image.  It  is  possible  to  “refocus”  the  SAR  image  so  that  all 
moving  targets  with  £=  £f  at  a  given  range  y=yF  appear  well  focused.  Refocusing  is 
accomplished  by  applying  a  filter  with  the  following  transfer  function; 


Jkl+epkl 


(54) 


Applying  the  filter  in  Equation  (54)  to  the  moving  target  spectrum  in  Equation  (45)  produces  the 
point  target  spectrum 
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(55) 


iF{K,ky\ep,yp)  =  U{k^,ky\£p,yp)  i{k^^,ky)  =  e 

provided  that  e=  Ef  and  ys-yp-  If  the  filter  parameters  Ef  and  yF  differ  from  the  target 
parameters  f  and  focus  quality  will  degrade. 

The  focusing  filter  may  be  applied  by  computing  the  two-dimensional  FFT  of  the  image 
data,  multiplying  by  the  transfer  fimction  Eqiiation  (54),  and  computing  the  inverse  FFT.  To 
provide  exact  focus  at  different  ranges  in  the  image,  the  entire  process  must  be  repeated  with 
different  values  ofyf- 


2.4.5  Refocusing  Moving  Targets  Via  Frequency-Domain  Resampling 

The  method  described  in  the  previous  section  will  correctly  focus  moving  targets,  but  may 
require  many  filtering  operations  to  cope  with  the  range-dependence  of  the  filter.  This  section 
will  describe  an  alternative  refocusing  method  that  provides  exact  focus  for  all  ranges 
simultaneously.  The  basic  idea  is  to  resample  die  frequency  domain  image  data  as  follows: 


■Spkl 


(56) 


The  coordinate  transformation  in  Equation  (56)  is  analogous  to  the  Stolt  transformation 
performed  in  the  Ctpk  image  formation  algorithm  [8].  Applying  Equation  (56)  to  the  moving 
target  spectrum  in  Equation  (45)  yields 


)^JC 


-“A  ■^<4+(e-ep)kl 


(57) 


When  =  £;  the  target  is  perfectly  focused: 

=  (58) 

There  is  one  important  subtlety  in  Equation  (56)  that  must  be  addressed  for  practical 
applications.  The  image  spectrum  I  in  Equation  (56)  is  the  Fourier  transform  of  the  entire 
(infinite)  image  plane.  In  practice  one  works  with  ,  the  transform  of  a  finite-sized  array  of 
image  data  centered  on  some  point  (xcyc)-  The  two  are  related  as  follows: 

i‘‘‘‘“‘{k^,ky)  =  i{k^,ky)e^‘=’"^*'^’='"^  .  (59) 


21 


The  same  relation  also  applies  to  the  refocused  image.  Combining  Equations  (56)  and  (59)  leads 
to  the  result 


rdata 


^yci^y 

o  ^ 


(60) 


iyj  ky-Jky-efk^ 

Comparing  Equations  (60)  to  (56)  shows  that  the  extra  phase  factor  e  ^  ^  is 

required  for  focusing  .  The  refocusing  process  in  Equation  (60)  becomes  especially  easy  to 

j  j  ‘^data  **data 

understand  when  rewritten  as  a  two  step  process  /  -4  /j  ->  /^r  ,  where 

=  (61) 

[k„ky',ep)  =  7/'"" {k„^k] -Spkl )  .  (62) 


The  first  step  in  Equation  (61)  applies  the  focusing  filter  U  specified  in  Equation  (54);  it 
provides  exact  focus  for  targets  at  range  y  =  y^.  The  second  step  in  Equation  (62)  is  a  non-linear 

resampling  in  the  ky  direction,  which  extends  exact  focus  to  all  range  values. 
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3.  (D-K  IMAGE  FORMATION 


This  section  will  demonstrate,  within  the  stationary  phase  approximation,  that  the  co-k 
image  formation  equation  produces  stationary  and  moving  target  signatures  precisely  identical  to 
those  produced  by  backprojection.  A  direct  calculation  for  the  stationary  target  case  will  be 
omitted,  since  it  is  just  a  special  case  of  the  moving  target  result. 

3.1  IMAGE  FORMATION  EQUATION 

The  basic  input  to  image  formation  is  the  range  compressed  data  D(^t).  The  co-k  equation 
computes  the  image  spatial  spectrum  directly: 


I{k^,ky)  =  \ 


OO  OO 


-ikj^-4yCkt 


for  Q, 


0 


otherwise. 


(63) 


The  band  support  region  Q  is  defined  according  to  Equation  (13).  As  in  previous  sections, 

k  =  yjk^  +ky  .  The  image  is  computed  by  an  inverse  Fourier  transform  of  I .  More  detailed 

discussion  of  the  (O-k  algorithm  (also  known  as  range  migration)  and  its  history  may  be  found  in 
[8].  In  the  published  literature,  there  are  some  variations  in  the  multiplicative  factors  included  in 
Equation  (63).  The  particular  factors  used  above  are  chosen  to  yield  an  ideal  bandpass  transfer 
fimction. 

3.2  STATIONARY  PHASE  ANALYSIS  (MOVING  TARGET) 

Consider  a  point  scatterer  target  moving  at  constant  velocity,  as  modeled  in  Section  1.2.3. 
The  range  compressed  data  is  strongly  peaked  at  t  =  lR^i^lc,  which  enables  the 

following  approximation: 


(64) 

A  more  rigorous  approach  to  the  approximation  in  Equation  (64)  is  to  regard  it  as  the  first  term  in 
a  Taylor  expansion  of  i^''^  around  the  point  t  =  2R^((^/c.  A  rather  laborious  calculation,  not 
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reproduced  here,  shows  that  the  correction  to  the  image  from  the  next  term  in  the  Taylor 
expansion  is  smaller  by  a  factor  of  order  IjkR  ~  XjR . 

Inserting  the  moving  target  data  into  Equation  (63)  and  applying  the  above  approximation 
yields 


hK^K)  = 


yjlTtk 


jdf  \dt  D,(.M 


(65) 


To  simplify  the  notation,  it  is  assumed  that  ikj^,ky)&Q..  The  integral  over  t  simply  produces  the 
Fourier  transform  of  the  range  compressed  data: 


\dl  c 


-^ckt  ^ 


(66) 


Making  the  change  of  variable  ^  ,  Equation  (65)  becomes 


where  the  phase  function  ^  is  defined  as 

0(0  =- 

■  -k,-kJ^ 

7 

The  derivatives  of  the  phase  function  are 

ao  k^  kc 

A 

and  - ; 

dO 

+7; 


kC 


(67) 


(68) 


(69) 


The  stationary  phase  condition  ||  =  0  is  satisfied  at  the  point  ^  ^ ,  where 


(70) 
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The  modified  band-support  region  is  defined  in  Equation  (39);  it  incorporates  the  condition 

^<y  required  for  the  existence  of  the  stationary  point.  As  promised,  the  image  spectrum 
in  Equation  (72)  is  identical  to  the  backprojection  spectrum  in  Equation  (45). 
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4.  ANTENNA  AND  TARGET  PATTERN  EFFECTS 


In  previous  sections,  it  was  assumed  that  the  SAR  antenna  has  a  transmission  pattern 
uniform  in  frequency  and  direction.  Similarly,  it  was  assumed  that  the  target  and  terrain 
scattering  patterns  are  uniform.  In  this  section,  the  uniformity  assumptions  will  be  relaxed, 
allowing  arbitrary  (but  slowly  varying)  patterns  instead. 

Figure  6  shows  the  basic  platform/target  geometry.  In  slant  plane  coordinates,  the  platform 
position  is  (^,0)  and  the  target  position  is  (xm,ym).  The  angle  between  the  flight  track  and  the  line- 
of-sight  to  the  target  is  at.  The  patterns  are  expressed  in  terms  of  the  complementary  angle 
d  =  7tll-K  and  the  frequency  f.  Note  that  0is  also  the  angle  with  respect  to  broadside.  The 
SAR  antenna  pattern  will  be  denoted  by  the  complex-valued  function  Aant(f,0),  and  the  target 
scattering  pattern  by  Atarg(f,d).  Both  patterns  also  depend  implicitly  on  the  cross-track  position 
of  the  target. 

Xm-^ 


<■ 


ym 


/ 

/ 

/ 


/ 

/ 


/ 

/ 

/ 

/ 


/ 

/ 


1/ 


Figure  6,  Platform/target  geometry. 


The  angle  ^is  related  to  the  platform  position  by  the  equation 


^  =  tan'' 


I  / 


(73) 
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The  basic  range  compressed  data  of  Equation  (3)  is  easily  modified  to  including  the  antenna  and 
target  pattern  factors: 


1  /max 
^  /min 

4.1  STATIONARY  TARGETS 


For  the  stationary  target  case, 


(^m»Tm)  =  (^0»To) 


and 


0  =  tan' 


^  To  ) 


(74) 


(75) 


The  backprojection  stationary  phase  analysis  in  Section  2.2  may  be  repeated  using  the  modified 
range  compressed  data  in  Equation  (74).  The  antenna  and  target  pattern  factors  propagate 
directly  through  the  calculation  to  the  final  result.  The  stationary  point  coordinates  are 
unchanged  since  the  phase  fimction  O  is  not  affected  by  the  slowly  varying  pattern  factors.  The 
stationary  point  values  for / and  6  are 


4;r  4;r 

The  backprojection  image  spectrum  is 


/  =—k=—Jk^  +kl  and  ^  =  tan"’—  . 

Arr  Arr  V  ^  K 


i(k  ,k  )  =  J  {or(k..k,)ea 

^  0  iov  (k^,k  Q, 


''X*  y  > 


where  the  antenna  and  target  transfer  functions  are  given  by 

Han,(K^ky)  =  A^„, 


^  c  k  ^ 

— tan"  — 
^4;r 


( 


^targ  >  ^>> )  ^targ 


.-1 K 


\ 


— A:,tan"‘^ 
^4;r 


(76) 


(77) 


(78) 


In  short,  the  antenna  and  target  patterns  are  manifested  as  filters  that  act  on  the  slant-plane  image 
data.  Note  that  the  relation  between  the  /Ts  and  the  yf’s  is  just  a  transformation  between  polar 
coordinates  and  Cartesian  coordinates. 
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Despite  the  identical  mathematical  properties  of  the  antenna  and  target  patterns,  there  is  an 
important  distinction  in  terms  of  practical  application.  The  target  pattern  will  differ  from  target 
to  target,  providing  useful  information  characterizing  the  targets.  The  antenna  pattern  is  fixed, 
and  it  obviously  provides  no  information  about  targets.  Any  nonlinear  behavior  in  the  phase  of 
(k^ ,  ky )  will  cause  blurring  in  the  image. 

To  generate  calibrated  imagery,  a  post-image  formation  filter  must  be  applied  to 
compensate  for  the  antenna  pattern  transfer  function.  The  calibration  filter  is  defined  by 
inverting  the  transfer  function  of  the  antenna: 


1 

^ant  ) 


(79) 


In  practice,  the  inversion  will  be  approximate  for  several  reasons.  First,  an  FIR  digital  filter 
needs  to  have  a  finite  number  of  taps.  Second,  the  exact  inversion  in  Equation  (79)  can  become 
singular  or  nearly  singular,  which  would  result  in  undesirably  large  amplification  of  the  noise 
floor.  Finally,  knowledge  of  the  antenna  pattern  is  almost  certainly  approximate,  so  it  is 
pointless  to  attempt  an  exact  inversion. 

The  calibrated  image  spectrum  is 


ici.K,ky) 


for(k„ky)eQ 
0  for(k^,ky)?:Q. 


4.2  MOVING  TARGETS 


(80) 


The  moving  target  coordinates  (x„,y„,)  are  given  in  Equation  (7).  The  backprojection 

stationary  phase  analysis  in  Section  2.4.1  may  be  repeated  using  the  modified  range  compressed 
data  in  Equation  (74).  The  antenna  and  target  pattern  factors  propagate  directly  through  the 
calculation  to  the  final  result.  The  stationary  point  coordinates  are  unchanged  since  the  phase 
function  <I>  is  not  affected  by  the  pattern  factors.  The  stationary  point  values  for /and  0  are 

f  ^ 

-1 


The  angle  y/  describes  the  direction  of  the  target  velocity  relative  to  the  radar  platform: 


f  =—k=—Jkl  +  kl  and  6-yr  +  tm 

4^  4^  ^ 


r^k^y+sk^ 


(81) 
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(82) 


*  -1 

yf  =  taxi  — 

\-u^ 

The  normalized  target  velocity  components  {u^,Uy)  can  be  expressed  in  terms  of  the  motion 
parameters  y  and  y/: 


\-u^  =ycosy/ 
Uy  =  j'sin  y/. 


(83) 


The  baclq>rojection  image  spectrum  (prior  to  antenna  pattern  calibration)  is 

iiK,ky)  =  AarXf^^)^an,Cf,G)i„opatiK,ky)  ,  (84) 

where  I„opaAK’^y)  is  the  moving  target  spectrum  without  pattern  effects,  as  given  in 

Equation  (45).  The  pattern  factors  may  be  written  in  terms  of  the  transfer  functions  defined  in 
Equation  (78): 


ant 


—cosw+J, 

.  y 


kl+ekl  sin^/, 


+skl  cosy/ — 


-sm^ 


(85) 


An  analogous  equation  hold  for  . 

Since  the  antenna  pattern  transfer  function  is  modified  by  target  motion,  calibration 
becomes  more  complicated.  In  general,  the  calibration  filter  in  Equation  (79)  will  not 
compensate  for  the  transfer  function  in  Equation  (85).  The  only  exception  is  when  the  antenna 
pattern  is  independent  of  the  direction  0.  Target  motion  only  affects  the  stationary  point  value 
of  not  /,  so  the  antenna  pattern  transfer  function  is  the  same  as  for  stationary  targets.  The 
calibration  filter  of  Equation  (79)  can  be  used  to  remove  the  antenna  pattern  effect,  and  then 
moving  targets  may  be  focused  as  described  in  Sections  2.4.4  and  2.4.5. 

4.3  MOVING  TARGET  REFOCUSING  AND  CALIBRATION 


The  moving  target  refocusing  algorithm  of  Section  2.4.5  may  be  applied  as  is  to 
uncalibrated  image  data  with  nonimiform  target  and  antenna  patterns.  In  other  words,  the  image 
data  given  by  Equation  (84)  can  be  plugged  directly  into  the  refocusing  formula  of 
Equation  (56).  For  a  point  target  located  at  (xo,^^) ,  the  result  is 
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(86) 


k  k 

—cosyf+k^smw,  k^cosw — ^sin^ 

KY  '  r  j 


-iXsk^-iy^ky 


where 


H(k,,k^)  =  H^,(k:c,ky)H^ik,,k^) 

is  the  lumped  transfer  function  of  the  antenna  and  target  patterns.  Note  that  Equation  (86) 
assumes  that  the  correct  value  of  £  was  used  when  applying  the  refocusing  equation. 

The  antenna  pattern  in  the  refocused  image  can  be  removed  by  applying  a  calibration  filter 
with  the  transfer  fimction^ 


H^ik„ky)  =  H 


/'  k  k  ^ 

—cosw+kySinw,  k^cosw — -sinw 

U  ''  7  J 


(87) 


A  key  observation  is  that  only  one  parameter  (f ,  or  equivalently,  j^)  is  required  to  focus  moving 
targets,  but  the  extra  parameter  y/  is  required  for  antenna  pattern  calibration. 

Even  after  applying  antenna  calibration  with  the  correct  value  of  yr.  Equation  (86)  shows 
that  the  target  pattern  is  still  distorted  by  a  mix  of  rotation  and  stretching  in  the  (k^,k^)  plane.  It 

is  possible  to  obtain  imdistorted  target  patterns  by  using  an  alternative  calibration  procedure. 
Rather  than  applying  the  filter  in  Equation  (87),  the  following  transformation  can  be  applied  to 
the  refocused  image: 


£,y/)  =  ip{Yk^cosy/-'ykySmyf,  k^smy/->rkyCo&y/\e).  (88) 


For  the  point  target  case,  the  result  is 


The  coordinates  (x^,  jyrp)  are  “fixed-time”  coordinates  of  die  moving  target.  In  other  words, 
they  are  the  target  coordinates  at  the  time  the  platform  reaches  the  position  ^  =  0  along  the 
synthetic  aperture.  They  are  related  to  the  broadside  coordinates^  (^o>^o)  ^  follows: 


^  The  inversion  is  to  be  understood  with  the  same  caveats  that  apply  in  the  stationary  target  case. 
*  i.e.,  the  position  of  the  target  when  it  is  broadside  to  the  SAR  platform;  see  Equation  (7). 
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Xjrj.  =  JToO  -  =  XoTCOS  y/ 

yFT=yo-  =^yQ-  y/. 


(90) 


The  transformed  image  data  in  Equation  (89)  can  then  be  calibrated  using  the  ordinary 
stationary-target  calibration  filter  in  Equation  (79).  The  target  pattern  is  undistorted  by  virtue  of 
the  transformation  in  Equation  (88).  As  a  side  benefit,  the  target  image  appears  at  the  true  fixed¬ 
time  location  {xpj,ypj) ,  instead  of  the  shifted  position  (x^,y,) . 

There  are  two  important  points  practical  points  that  are  worth  addressing  here.  The  first 
point  concerns  the  fact  that  true  image  data  is  typically  centered  on  some  point  (x^jj^)?  while 

the  above  equations  implicitly  assume  the  image  data  is  centered  at  zero.  The  matter  was 
addressed  in  Section  2.4.5,  and  the  same  solution  applies  here.  The  second  point  is  that  the  two- 
dimensional  resampling  step  in  Equation  (88)  would  be  easier  and  more  efficient  to  implement  in 
terms  of  two  sequential  one-dimensional  resampling  operations. 

Taking  the  above-mentioned  points  into  account,  a  multi-step  procedure 

rdata  ^data  ^data  -^data  ’••data 

I  is  suggested  for  moving  target  focusing  and 

antenna  calibration.  The  steps  are  as  follows: 


(  k 

[k^,ky',  ep,y/ir)  =  if“‘  k^,-^Xmy/p+kyStcy/p-,  Sp 

i‘^‘‘[k^,ky-,  £p,y/p)  =  I'^“‘(yk^cosy/p-YkyS\nyfp,  k^-,  £p,y/p) 

^ant 


(91) 

(92) 

(93) 

(94) 

(95) 


The  “data”  superscript  has  the  same  meaning  as  in  Section  2.4.5,  indicating  image  data  centered 
at  {x^,yc)  •  The  first  two  steps  in  Equations  (91)  and  (92)  are  precisely  identical  to  the  focusing 

algorithm  given  in  Equations  (61)  and  (62).  The  next  two  steps  in  Equations  (93)  and  (94) 
perform  one-dimensional  resampling  operations  that  undo  the  scaling  and  rotation  effects  in  the 
antenna  and  target  patterns.  Note  that  the  combined  effect  of  Equations  (93)  and  (94)  is  identical 
to  that  produced  by  Equation  (88),  but  the  latter  requires  a  more  complicated  two-dimensional 
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resampling.  Finally,  Equation  (95)  calibrates  the  image  by  eliminating  the  antenna  transfer 
function.  Caveats  discussed  in  the  previous  section  regarding  inversion  of  the  antenna  transfer 
function  apply  equally  well  here. 

Applying  Equations  (91-95)  to  a  point  target  with  the  correct  values  £f=£  and 
produces  the  spectrum 


jdata 


^kx,ky,  B^y/^  —  Hi^{kxiky)6 


-ixsFTkx-VSFTky 


The  target  is  focused  at  the  “shifted  fixed-time”  coordinates  (x^py,  j^pj.)  that  are  related  to  the 
“fixed-time”  coordinates  in  Equation  (90)  as  follows: 

XSFT  =  ^FT  +  x,(l-  r  COS  i/f)  -y^smi/F 

ysFT  +x^2'sun/+y^(l— cos^). 

A  nice  feature  of  the  multistep  focusing  procedure  above  is  that  the  operations  occur  in 
order  of  importance  (most  to  least)  as  well  as  order  of  complexity  (least  to  most).  Applying  the 
first  step  in  Equation  (91)  requires  just  one  focusing  parameter,  no  resampling,  and  will  mitigate 
most  of  the  blurring  due  to  target  motion.  The  second  step  introduces  no  new  parameters, 
requires  only  one-dimensional  resampling,  and  extends  focus  to  all  ranges.  The  remaining  steps 
require  a  second  parameter  and  multiple  resampling  steps  to  achieve  an  undistorted  target 
response  with  the  antenna  pattern  removed.  The  number  of  steps  actually  implemented  in  a 
practical  system  can  be  chosen  to  suit  the  desired  accuracy  and  available  computing  resources. 
Conceptually,  each  additional  step  makes  a  correction  at  a  more  subtle  level  than  the  previous 
one,  with  the  full  sequence  providing  an  exactly  focused  and  calibrated  image. 
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APPENDIX 

STATIONARY  PHASE  EQUATIONS 
A.1  ONE  DIMENSIONAL  INTEGRALS 
Consider  an  integral  of  the  form 


l=\g{x)e'f^^^^dx,  (A.1) 

V 

where  ^x)  is  strictly  real-valued.  As  //  becomes  large,  the  phase  factor  will  fluctuate 
rapidly.  The  rapid  fluctuations  tend  to  cancel  out,  except  at  stationary  points  where  -^  =  0. 
Asymptotic  analysis  in  the  limit  //  -^  +«>  leads  to  the  stationary  phase  equation 


lit 


(A.2) 


The  sum  is  over  the  set  of  stationary  points  x^&V  where  <l>\x^)  =  ^\x=xk 

stationary  points  in  the  integration  region  V,  then  / « 0  at  leading  order  in  the  asymptotic 
expansion.  The  “sgn”  function  in  Equation  (A.2)  simply  returns  the  sign  of  its  argument: 


f+l  ifx>0 


sgn(x)  =  «^ 


0 

-1 


ifx  =  0 
ifx<0. 


(A.3) 


Equation  (A.2)  applies  under  “generic”  conditions:  ^  0  and  no  stationary  points  on 

the  boundaries  of  the  integration  region  V.  A  careful  derivation  of  (A.2)  and  discussion  of 
conditions  on  the  functions  and  g(x)  may  be  found  in  a  variety  of  references;  an  excellent 
exposition  is  provided  in  [9].  For  the  purposes  of  this  document,  the  integrals  evaluated  are 
always  of  “generic”  type  and  Equation  (A.2)  will  be  applied  without  further  discussion  of  any 
subtleties. 

The  derivation  of  Equation  (A.2)  may  be  sketched  as  follows.  At  leading  order,  the  only 
parts  of  the  integral  Equation  (A.1)  that  matter  are  near  the  stationary  points.  In  odier  words, 

4+f  s 

/  «  S  J  =  5^  J g(4  +  (A.4) 

k  x/i-s  k  -e 
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for  some  suitable  value  of  £  Expanding  (f>  to  second  order  aroimd  the  stationary  point(s)  and 
treating  g  as  a  constant  yields 


I  =  Sirfty***'  & .  (A.5) 

k 

Since  the  integrand  fluctuates  rapidly  for  |5|  >£,  the  limits  of  integration  may  be  extended  to 

infinity  without  changing  results  at  leading  order.  The  resulting  Gaussian  integral  may  be 
computed  exactly,  producing  the  final  result  in  Equation  (A.2). 

A.2  MULTIDIMENSIONAL  INTEGRALS 

The  iV-dimensional  stationary  phase  integral  is  of  the  form 

(A.6) 

V 

The  bold  S5mibol  x  denotes  an  N-dimensional  real  vector,  and  denotes  integration  over  the  N 
components  of  x.  The  stationary  phase  derivation  proceeds  exactly  as  in  the  one-dimensional 
case.  The  stationary  points  are  x^^  e  F  such  that  VO(x^ )  =  0 . 

After  expanding  to  second  order  around  the  stationary  pomt(s),  the  integral  becomes 

/  «  .  (A.7) 

The  notation  represents  the  matrix  of  second  derivatives  of  <I>  at  the  k-th  stationary 

point.  Performing  the  Gaussian  integral  in  Equation  (A.7)  leads  to  the  TV-dimensional  stationary 
phase  equation: 


7 

*  deta^4.(ii)  ' 


(A.8) 


The  matrix  signature  function  in  Equation  (A.8)  is  a  matrix  generalization  of  the  “sgn”  function: 


sig(M)  =  (#  of  positive  eigenvalues  of  M)— (#  of  negative  eigenvalues  of  M) . 
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